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Diffusion is a key element of a large set of phenomena occurring on natural and social systems 
modeled in terms of complex weighted networks. Here, we introduce a general formalism that allows 
to easily write down mean-field equations for any diffusive dynamics on weighted networks. We also 
propose the concept of annealed weighted networks, in which such equations become exact. We 
show the validity of our approach addressing the problem of the random walk process, pointing out 
a strong departure of the behavior observed in quenched real scale- free networks from the mean-field 
predictions. Additionally, we show how to employ our formalism for more complex dynamics. Our 
work sheds light on mean-field theory on weighted networks and on its range of validity, and warns 
about the reliability of mean-field results for complex dynamics. 
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I. INTRODUCTION 

Weighted networks represent the natural framework 
to describe natural, social, and technological systems, in 
which the intensity of a relation or the traffic between el- 
ements are important parameters In general terms, 
weighted networks (WN)are an extension of the concept 
of network or graph 0, 3 , in which each edge between 
vertices i and j has associated a variable, called weight. 



taking the form 



g{i,j)aij, where is the ad- 



jacency matrix and g{i,j) is a function of i and j 
Practical realizations of weights in real networks range 
from the number of passengers traveling yearly between 
two airports in the airport network , to the intensity of 
predator-prey interactions in ecosystems or the traffic 
measured in packets per unit time between routers in the 
Internet Q . A non- weighted network can thus be under- 



stood as a binary network, in which 



taking 



the value 1 when vertices i and j are connected, and zero 
otherwise. Leaving aside the issue of the proper char- 
acterization of the topological obscrvables associated to 
networks and their weighted extensions P, 0, [^, l9l-[Tl| , a 
most relevant aspect of WN is the effect of the weight 
distribution and structure on the properties of dynami- 
cal processes taking place on top of them [l^ - [l9| . These 
can have an impact in problems ranging from the in- 
formation transport in mobile communication networks 
[20j . the risks of congestion on Internet 21 1, the disease 
spreading in the air transportation network [22j or bi- 
ological issues such as the role of weak interactions in 
ecological networks [23j . 

In general, the theoretical understanding of dynamical 
processes on networks is based at a mesoscopic level on 
the heterogeneous mean-field (HMF) app roximation and 
the annealed network approach |24l . |25| , which assume 
that the degree is the only relevant variable character- 
izing the vertices, and that dynamical fluctuations arc 
not relevant. All dynamical variables can therefore be 
described in terms of deterministic rate equations, as a 
function of time and degree. Here, we develop the HMF 
theory for general local diffusive processes on WN by ex- 



tending the annealed network approximation [2J] to the 
weighted case. We also deflne the concept of annealed 
weighted network, in which HMF theory is exact. We 
show the validity of our approach by obtaining an exact 
mean-field solution of the paradigmatic random walk pro- 
cess on WN. A comparison with numerical simulations 
allows to set the limits of validity of HMF on quenched 
weighted networks, showing also that it yields in some 
cases severely incorrect results. Finally, we present an 
example of the application of the annealed weighted net- 
work approximation to more complex dynamics on WN 
to obtain HMF rate equations. 

Our approach allows for a straightforward mean-field 
description of diffusive dynamics on weighted networks, 
but at the same time calls for a careful interpretation 
of mean-field results on quenched networks, an issue al- 
ready discussed in the non- weighted case j26l - l29j . HMF 
indeed does not take into account the quenched structure 
of the considered network, which however can play a fun- 
damental role inducing strong deviations even in the ap- 
parently safe case in which only degree dependent weights 
are present. This should be recalled when mean- field pre- 
dictions are to be compared with numerical simulations 
or used as a baseline to make predictions on processes on 
real weighted networks. 



II. ANNEALED WEIGHTED NETWORK 
APPROXIMATION 

A. Annealed network approximation in binary 
networks 

The HMF theory for dynamical processes on binary 
complex networks is ultimately based on the annealed 
network approach [2^ . that consists in replacing the ad- 
jacency matrix a^ , which in a real quenched network is 
composed by zeros and ones at fixed positions, by the 
probability aij that vertices i and j are connected. At 
a statistical degree coarse-grained level, we can replace 



by the ensemble average a{ki, kj) [30[, defined as the 



2 



probability that two vertices of degree ki and kj are con- 
nected; that is, by the degree class average (for undi- 
rected networks) 



a{k,k') 



1 1 _ kP{k'\k) 



NP{k) NP{k'\ 



NP{k') 



where i g fc denotes a sum over the set of vertices of de- 
gree k, P[k) is the degree distribution (probability that 
a randomly chosen vertex is connected with other k ver- 
tices) 3 and P{k'\k) is the conditional probability that 
a vertex of degree k is connected to a vertex of degree k' 
[3ll |. In other words, the original network is replaced by a 
fully connected graph in which each edge has associated 
a connection probability a{ki,kj) that depends only on 
the degree of its endpoints. The annealed network ap- 
proach can be directly applied in calculations at the mi- 
croscopic (vertex) level or, otherwise, fruitfully extended 
to a more intuitive mesoscopic (degree class) level. In the 
case of diffusive dynamics, the reasoning runs as follows. 
In diffusive dynamics, particles or interactions move from 
one vertex to another by jumping to a randomly chosen 
nearest neighbor. Therefore, the probability that vertex 
i interacts with vertex j, P{i — > j) (the propagator of the 
diffusive interaction), is given by 



this case, the annealed network is completely defined by 
the quantities P(k) and P{k'\k). From a computational 
point of view, a diffusive dynamics is easily simulated in 
these networks, by assigning to each vertex of degree k 
nearest neighbors that are randomly selected among all 
the vertices in the network, according to the probability 
P{k — >■ k') = P{k'\k). In the case of degree uncorrelated 
networks, this interaction probability simplifies further 
to k'P{k')/ (fc), so that nearest neighbors are randomly 
chosen with probability proportional to their degree. Due 
to the very definition of the annealed network approach, 
HMF theory is exact in annealed networks [soj . 

B. Extension to weighted networks 

We now generalize the annealed network approach to 
the weighted case. Assuming that the diffusion process 
is local (depends only on the departing and arriving ver- 
tices), the most general way to take into account weights 
in the random walk is to define the probability that a 
vertex i interacts with vertex j as a general function of 
the weight Wy. The normalized form of this probability 
will be 



(4) 



P(z ^ j) = 



k. 



(2) 



Performing a coarse-graining in degree classes, we can 
define the probability that a vertex of degree k interacts 
with a vertex of degree fc', namely 



P{k k') 



Pik'\k), 



(3) 

where we have used Eq. (H]). With this concept in hand, 
writing down the HMF rate equations for any diffusive 
dynamical process turns out to be straightforward [23 |. 

Equation Q also encodes the idea of an annealed net- 
work [29, 30, 32 In opposition to quenched net- 
works, in which edges are frozen between pairs of vertices, 
annealed networks are dynamical objects, changing in 
time over a certain timescale tjv, while keeping constant 
the degree distribution P{k) and correlations P{k'\k). 
When the network time scale is very small, tjv — > 0, 
connections in the network are completely reshuffied be- 
tween any two microscopic steps of the dynamics. In 



The only restriction we impose to the function F[x] is 
that F[0] = 0, in order to avoid diffusion between dis- 
connected vertices (i.e., with aij ~ 0). This implies that 
we can write 



F[aj,,] = F[g{i,j)]a., 



(5) 



At a mesoscopic level, we again consider the average 
probability of interaction of a vertex k with a vertex fc', 
that is defined, in analogy with Eq. by 

[7vp(fc)l-iy;.^, v.^,, F\uj,A 

[NP{k)] J2iek T,r PM 

At the HMF level, where the vertices's degrees are the 
only relevant topological variables, this expression can be 
simplified considering that the dependence of the weight 
on the vertices at the endpoints of each edge can be ex- 
pressed as a function of the corresponding degrees, that 
is, LUij = g{ki, kj)aij . Using this relation, we have from 
Eqs. © and (P 



3 



From this expressions, wc finally obtain the degree 
coarse-grained weighted propagator 



F[g{k,k')]P{k'\k) 
E,F[5(fc,q)]P(q|fc)' 



(9) 



satisfying the normalization condition ^j,, Pio{k — > k') ~ 
1. 

Simpler expressions for the weighted propagator can be 
obtained if we consider a linear diffusion process [H, , 
proportional to the weight w^-, in which F{x) = x. Fur- 
thermore, considering weights that are symmetric multi- 
plicative functions of the degrees at the edges' endpoints, 
namely gik-, k') = gs{k)gs{k') (as is the case, for example, 
in the airport transportation network 0), and a degree 
uncorrelated network, with P{k'\k) = k'P(k')/ (fc), we 
are led to the simple expression 



Pn,{k^k') 



k'g,{k')P{k') 
{kUk)) 



(10) 



where we define (G(fc)) = J2k G{k)P{k). 

Writing the mean-field rate equations for most diffusive 
dynamical processes on WN is now easy: The analytics 
describing the dynamics on simple binary networks can 
be generalized to the weighted case simply by replacing 
the propagator P{k k') with Pw{k — )• k'). 

The spirit in which the HMF rate equations are con- 
structed leads also to introduce the concept of an an- 
nealed weighted network. In this case, in the event of an 
interaction, a vertex of degree k chooses as the interaction 
target a vertex of degree k' , randomly selected among 
all vertices in the network with probability Pw{k — > fc'). 
Therefore, an annealed weighted network can be under- 
stood as an annealed binary network, which is completely 
reshuffied between dynamic time steps, preserving its de- 
gree distribution P{k) and an effective degree correla- 
tion pattern given by Pes{k'\k) = P^{k k'). Again, 
in the case of linear diffusion in a degree uncorrelated 
network with symmetric multiplicative weights, this con- 
nection probability simplifies to k'gs{k')P{k')/{kgs{k)), 
so that the interacting vertex is chosen with probability 
proportional to kg^ (fc) . Simulations run on such networks 
are described exactly by weighted HMF rate equations. 
Quenched WN can in principle determine different be- 
haviors, and, as we will see below, this is in fact the case. 



III. RANDOM WALKS ON WEIGHTED 
NETWORKS 

A. Weighted heterogeneous mean-field solution 

To check the application of the formalism derived 
above, let us consider as an example the simple yet 
paradigmatic random walk process |35| . This is defined 
by a walker that, located on a given vertex of degree k at 
time t, hops to one of the k neighbors of that vertex at 
time t + 1, randomly chosen with probability depending 



on the weights connecting the vertices. The simplicity of 
this problem allows to solve it exactly for certain forms of 
the weight structure [13, , following the master equa- 
tion approach developed in Ref. [35| . Here we focus on 
the application of the simpler HFM in the most general 
case, which provides results for any weight pattern. 

To construct the appropriate HMF rate equation, we 
consider the probability Rk{t) that the walker is in any 
vertex of degree fc, which fulfills the master equation 

Rkit) = -J2Rk{t)P^{k^k')+Y,PUk'^k)Rk'{t) 

k' k' 

= Rk{t)+Y,P^ik' ^k)Rk'{t). (11) 

k' 

The first term in Eq. (|lip represents the outflow of prob- 
ability due to walkers abandoning vertices of degree fc, 
while the second term represents walkers arriving to ver- 
tices of degree k from vertices of degree fc', following a 
weighted random step. In the steady state this equation 
takes the iterative form 



Rk = J2PUk' ^k)Rk 



(12) 



Interpreted in terms of a Markov process 137| , the steady 
state solution is given, for any degree correlation and 
weight pattern, by the expression 



lim P,(,")(fc^fc), 



(13) 



where Pi,"'' is the n-th power of the probability Pw{k — >■ 
fc'), considered as a matrix. Explicit solutions in the 
steady state can be found for some weight patterns by 
imposing the detailed balance condition RkPwik — > k') = 
Rk'Pw{k' — >■ k) [13], which ensures a solution of Eq. (ffel) . 
From this condition we obtain, for a general diffusion 
process, 

Rk_ ^ Pn,{k' ^k) 

Rk' ~ PUk^k') 

^ F[g{k', k)]kP{k) Eg F[9i.k, q)]P{q\k) 

- F[g{k, k')]k'P{k') E, F[g{k', q)]P{q\k') 



,(14) 



where in the last equality we have used the degree de- 
tailed balance condition k' P{k')P{k\k') = kP{k)P{k'\k) 



[38l |. Then, for symmetric weights uji 



degree level, g{k,k') = gs{k,k') 
the normalized probability 



or, at the 



fRk 



gs{k',k), we obtain 
kP{k)j:^FMk,q)]P{q\k) 



Ek' E,k'P{k')F[Uk',q)]P{q\k')- 



(15) 



If the weights only depend on the first vertex in the edge, 
that is, uJij = gi{i)aij or g{k,k') = gi{k), the effect of 
weights effectively vanishes and we recover diffusion in a 
binary network |35|, 



Rk 



kP{k) 



Ek' k'P{k') 



(16) 
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Finally, if the weights only depend on the second vertex 
in the edge, that is, — g2{j)0'ij or g{k,k') = g2{k'), 
then we obtain 



Rk 



kP{k)FMk)]j:,FMq)]P{q\k) 

Ek' E„ k'Pik')F[-g2ik')]FMq)]PiqW) ' 



(17) 



Noticeably, results and ((T7| coincide if diffusion is 
linear, F{x) = x, and the weights are symmetric and mul- 
tiplicative, Q-tik . k') = gs{k)gs{k'), recovering the results 
reported in [ij, [3g . 

Further quantities can be analogously computed with 
ease in the annealed WN approximation. For exam- 
ple, the coverage S{t) [H, [H, H^], defined as the aver- 
age number of different vertices visited by the walker 



at time t, can be computed as follows. Let us de- 
fine Sk{t) = Vk{t) /[NP{k)] as the fraction of vertices 
of degree k visited by the random walker, Vk{t) be- 
ing total number of such vertices. From here, we have 
S{t) = NJ2k-P(k)sk{t). The quantity Vk{t) increases in 
time as the random walk arrives to vertices that have 
never been visited. Therefore, at a mean-field level, it 
fulfills the rate equation 

Vk{t) = [1 - Sk{t)]J2 Proik' ^ k)Rk'it). (18) 



Dividing this equation by NP(k) and using Eq. ([9]), we 
are led in the general case to 



hit) = [i-skmkY. 



F[g{k'MP{k'\k) 
Nk'P{k')j:F[-g{k',q)]P{q\k') ^ >■ 



(19) 



r 



Assuming now that Rk{t) reaches its steady state in a 
very short time, of order 1, we can substitute the steady 
state value Rk into Eq. (jl9|) and integrate it, with the ini- 



J 



tial condition s/c(0) — 0. For symmetric weights, Eq. (fT5)) 
the integrations leads to the result 



Sk{t) = 1 - exp I -t 



kY.k'F[9s{k'MP{k'W) 
NEk" E« k"P{k")F[g,ik",q)]P{q\k") 



(20) 



r 



Moreover, straightforward mean-field arguments |40l - |42| 
provide an expression for the mean first-passage time 
(MFPT) Tfc, defined as the average time that a walker 
takes to arrive to a given vertex of degree fc, starting 
from a randomly chosen initial vertex [35| . Let us define 
the occupation probability Pk{t) = Rk{t)/[NP{k)] as the 
probability that the walker is at a given vertex of degree 
k. Thus, the probability for the walker to arrive at a 
vertex i, in a hop following a randomly chosen weighted 
edge, is given by q(i) = Pkiit). Therefore, the probabil- 
ity of arriving at vertex i for the first time after t hops 
is Pa{i; t) = [l- q{i)Y-^q{i). The MFPT to vertex i can 
thus be estimated as the average 



(21) 



We have thus that the MFPT is proportional to the in- 
verse of the steady state probability, Rk, namely 



Tk 



NPik) 

Rk 



(22) 



This simple HMF approximation can be improved by the 
inclusion of sub-leading terms using more formal tech- 
niques [ssIIsgI. 



B. Dynamics on quenched weighted networks 

The weighted HMF expressions obtained above will ex- 
actly match, by definition, the results of simulations on 
annealed WN. A different issue is, however, their valid- 
ity for real quenched WN, in particular for networks with 
a scale-free degree distribution, P{k) ~ fc-^ [I,!!]. The 
question of the non mean field behavior of dynamical pro- 
cesses on binary networks has been already pointed out 
in the literature [1^ H^. As we will see, the situation 
can get even worse in the case of WN. 

We consider as a particular example the case of liner 
diffusion on symmetric multiplicative weighted networks, 
with weight intensity gs{k) = k^ , typical of some real sys- 
tems such as the airport network [6| . In the limit of large 
positive 6, in an annealed WN we would expect the ran- 
dom walker to reach the largest hub in the first time step. 
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FIG. 1: (Color Online) Top: Relative error between the the- 
oretical HMF prediction Eq. (|23p and numerical simulations 
for the MFPT in quenched scale-free WN of size N = lO" 
and 7 = 2.5. Bottom: Relative error for 9 = 2.5 and different 
network sizes. 



FIG. 2: (Color Online) Top: Average error as a function of 
9 in quenched WN of size N = 10* and 7 = 2.5 for the 
simple HMF estimate (circles), Eq. (|23|l . and a more refined 
approximation (crosses), Eq. (|27|l . Bottom: Range of 9 values 
yielding an average error E smaller than a given value for the 
HMF estimate. 



and then stay there forever [5l|. In a quenched WN, on 
the other hand, the walker is expected to commute for- 
ever between the first encountered pair of vertices who 
have the property of being reciprocally the highest degree 
neighbor of each other. In order to explore quantitatively 
this intuitive disagreement, we have performed numeri- 
cal simulations of random walks in degree uncorrelated 
and multiplicative scalc-frec WN, with weight intensity 
gs{k) — and linear diffusion, which yield the HMF 
results 



1 _ iV(fci- 



1 



TV 



(23) 
(24) 



where J-'-y^e{x) is a scaling function depending on 7 and 9. 
Uncorrelated scale-free quenched networks with any de- 
gree exponent 7 are created using the Uncorrelated Con- 
figuration Model (UCM) characterized by a hard 
cut-off kc = a/ZV, preventing the generation of correla- 
tions for 7 < 3. Random walks are performed by moving 
the walker from its present position i to a nearest neigh- 
bor j, chosen with probability 



kj 



J2q ^q^iq 



(25) 



On annealed WN, on the other hand, the next step of the 
walk is a vertex j, randomly chosen among all vertices 
in the network with probability proportional to ^j"*"^, see 
Eq. p^ . In our simulations we keep fixed a degree ex- 
ponent 7 = 2.5 and a minimum degree of the network 
TO = 4. 

We first study the validity of the weighted HMF predic- 
tion for the MFPT, by examining the differences between 



the theoretical expression r^'*, Eq. ((23)) . and the numer- 
ical result, r^"™. In particular, we focus on the relative 
error as a function of the degree, defined as 



E, = 



(26) 



Fig. HJtop) shows the numerical evaluation of this func- 
tion in quenched scale-free WN of fixed size N = 10^, 
for different values of 9 (H^]. For 9 = 0, corresponding 
to a binary network, we recover a constant baseline er- 
ror of about 15%, as already reported in the literature 
[4^ . As 9 increases, the numerical results start to de- 
viate from the HMF prediction, the error being larger 
for small degrees and also increasing for large k. The 
error in the estimate of the weighted HMF theory also 
increases when increasing the system size, as can be seen 
from Fig. [IJbottom). 

In order to provide an estimate of the range of va- 
lidity of the HMF approximation, in Fig. [5] we explore 
the average error as a function of 9, E = J2k 
As we can see from Fig. [Utop) , for the particular value 
7 = 2.5 considered, the error is minimal, compatible with 
the baseline of 15%, in the interval 9 € [0, 0.5], while it in- 
creases for values of 6* outside this interval. In Fig.[2Jtop) 
we additionally compare the performance of the HMF 
prediction with a more elaborate estimate of the MFPT 
obtained by the master equation technique, namely [11] 



N{k^+'>) N{ky^ 



k{k^ 



- 2. 



(27) 



As we can see, Eq. (P7)) fares slightly better for < 0, 
but the HMF result remains a fairly good first order ap- 
proximation. 

This kind of plot can be used to estimate the range of 
values of 9 for which the weighted HMF approximation 
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FIG. 3: (Color Online) Data collapse of the network coverage, 
as given by Eq. (|24[). in annealed (A) and quenched (Q) scale- 
free WN, with 7 = 2.5 for different values of Q and network 
sizes. Curves shifted on the horizontal axis for clarity. Inset; 
Unshifted curves for Q = 2.5 and A'^ = 10"' for annealed (^4, 
continuous line) and quenched (Q, dotted line) WN. 



provides a correct result within an accepted maximum 
tolerable error, see Fig. [2jbottom). Thus, for example, a 
deviation smaller than a 50% in scale-free net"works with 
7 = 2.5 is achieved for values of ^ € [—0.9, 1.2]. Outside 
this range, other more sophisticated approaches should 
be followed. 

Concerning the results for the network coverage <S'(t), 
a direct comparison with Eq. ([M|) is not feasible, since 
the exact form of the function J-^fi{x) depends on several 
approximations (steady state approximation in Eq. (|20|) , 
continuous degree approximations, etc.). We therefore 
take as the main weighted HMF prediction the scaling 
form in Eq. (j24[) , which we check by means of a data col- 
lapse analysis: If Eq. is correct, we expect that plots 
of 1 - S{t)lN as a function of t/{N{k^+'')) wiU collapse 
onto the universal function !F^^e{x) when plotted for dif- 
ferent values of N. In Fig. |3] we compare the results of 
simulations on quenched weighted networks with simu- 
lations in the corresponding annealed ones, for different 
values of 9. For small values of 9 (i.e. smaller than 0.5), 
we observe a perfect agreement in both sets of simulations 
with the weighted HMF prediction. Again, however, the 
perfect collapse of the annealed case contrasts with the 
large deviations shown by quenched networks for large 
values of 9, outside the regime of validity of the weighted 
annealed network approximation. 

The reason for the failure of weighted HMF theory in 
the random walk case is easy to understand. For positive 
9, a quenched topology reduces the accessibility of the 
network by trapping the walker in pairs of adjacent high 



degree vertices nodes and, at the same time, slows down 
the exploration of the network due to the trapping ef- 
fect of high degree vertices H^. These effects, which are 
stronger for large 0, explain the large MFPT observed for 
small and large degree. Fig. [Ijtop) , as well as the larger 
times needed to reach a fixed coverage. Fig. [31 For neg- 
ative 0, on the other hand, the walker is biased towards 
small degree vertices, missing thus the hubs that provide 
connectivity to far away regions in the network [46j. 



IV. GENERAL DYNAMICAL PROCESSES 

As we have pointed out in Sec. [Hi the annealed 
weighted network approximation can be easily applied to 
general diffusive processes in weighted networks to obtain 
an approximate HMF solution. Within this framework, 
the weighted HMF rate equations are obtained from the 
ones on binary networks by substituting the binary prop- 
agator by the weighted one Pw{k — J> k'), Eq. (|9]). 

As a case example of application, we will consider here 
the case of the contact process (CP) [4?! , w hose dynam- 
ics on a network is defined as follows [2g: An initial 
fraction pQ of vertices is randomly chosen and occupied 
by a particle. The time evolution of the process runs as 
follows: At each time step a particle in a vertex i is 
chosen at random. With probability p, this particle dis- 
appears. With probability 1 — p, on the other hand, the 
particle may generate an offspring. To do so, a vertex j, 
nearest neighbor of the vertex i, is chosen. If vertex j is 
occupied by a particle, nothing happens; if it is empty, a 
new particle is created on j. In any case, time is updated 
as, t ^ t + l/n(i), where n{t) is the number of particles 
present at the beginning of the time step. In a binary 
network, the second vertex j is chosen uniformly at ran- 
dom among all neighbors of i, i.e. with probability l/fc^. 
In a weighted network, on the other hand, it is selected 
with probability depending on ujij. 

The key point in the HMF analysis of the CP is the 
probability (t) that a vertex of degree k is occupied by a 
particle at time t [2^, [2^. To write the HMF equation for 
this quantity it is useful to focus instead on the number 
of occupied vertices of degree fc, rifc(t), which is related 
to pk{t) by 



Pk{t) 



nk{t) 
NP{k)' 



(28) 



The rate equation for nk(t) is easy to obtain by means of 
mean- field arguments |26l . l48j , leading to 



rifc (t ) = -puk {t) + {l-p)[l-pk{t)]Y^P^{k' 



k)nk' (t). 
(29) 

The first term comes from particles disappearing from 
vertices of degree k with probability p; the second term 
correspond to the offsprings of particles in connected ver- 
tices of degree k' (with probability 1 — p), arriving to 
vertices of degree k which arc empty (with probability 
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1-Pk{t)). Dividing Eq. (gH) by NP{k), substituting the 
expression of the weighted propagator Eq. ^ and apply- 



In this last equation, we have performed a rescaling of 
time, and defined the parameter A = (1 — p)/p. For 
symmetric multiplicative weights of the form g{k, k') = 
gs{k)gs{k'), and assuming linear diffusion a degree un- 
correlated network, we are led to the simplified equation 

puit) = -pu{i) + A[l - p,(t)]^||M, (31) 

where p{t) = P{k)pk(t). In this expression we recover 
the result presented in Ref. [T5| . 

Along the same lines, more complex dynamical pro- 
cesses subject to general weighted diffusive rules can also 
be considered and the corresponding weighted HMF the- 
ories easily developed. 

V. CONCLUSIONS 

In this paper we have developed a general formalism to 
write down HMF equations for diffusive dynamical pro- 
cesses on WN. Moreover, we have introduced the concept 
of annealed weighted networks, for which HMF theory 
represents an exact description. Considering as a simple 
example the random walk process, we have presented ex- 
act mean-field implicit solutions for its behavior for any 
correlation and weight pattern, as well as explicit for- 
mulas for the case of symmetric weights. By means of 
numerical simulations, we have also shown that weighted 
HMF theory can describe diffusion in real quenched scale- 
free networks with multiplicative weights proportional to 
k^ only for small values of \9\. Finally, we have demon- 
strated how the annealed weighted approximation can be 
fruitfully applied to obtain information for more complex 
dynamical systems, as for example the contact process, in 
which more exact analytical alternatives can be difficult 
(or impossible) to work out. 

Overall, our work provides a straightforward method 
to describe any dynamical process on weighted network 
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